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^ Abstract 

We consider the Bohr correspondence limit of the Schrodinger wave 
function for an atomic elliptic state. We analyse this limit in the con- 
^ text of Nelson's stochastic mechanics, exposing an underlying deter- 

, ^^ ^ ministic dynamical system in which trajectories converge to Keplerian 

motion on an ellipse. This solves the long standing problem of obtain- 
^^J ing Kepler's laws of planetary motion in a quantum mechanical setting. 

r — In this quantum mechanical setting, local mild instabilities occur in 

lO the Keplerian orbit for eccentricities greater than l/\/2 which do not 

— »^ occur classically. 

i> 1 Background 

o 

^ 1.1 Introduction 



One of the longest standing problems in quantum mechanics is to obtain Ke- 
pler's laws of planetary motion by taking the Bohr correspondence limit of 
a suitable Schrodinger wave function for an atomic elliptic state. The Bohr 
correspondence limit of a quantum state is found by allowing the quantum 
numbers to become infinite while h becomes vanishingly small with the ener- 
gies and angular momenta of the system fixed with physical values of order 
one [4]. We have approached this quantum mechanical problem from the 
perspective of Nelson's stochastic mechanics. 



In 1966, Edward Nelson produced a new formulation of non relativistic 
quantum mechanics in terms of diffusion processes satisfying the Nelson- 
Newton law, a stochastic version of Newton's second law of motion [12, 13, 
14] . This work was based on the assumption that every particle of mass m 
moves in a random environment driven by a Brownian motion with diffusion 
constant H/2m. Some mathematical physicists felt that Nelson's formulation 
merely recapitulated the standard results of Schrodinger quantum mechanics 
without explaining the physical agency responsible for the Brownian noise in 
the system, and so the theory did not prove popular. However, the theory 
does give a new mental picture of the quantum world in which one can calcu- 
late the statistics of new observables, like first hitting times, which should be 
open to measurement [8]. Moreover, within this formulation. Nelson showed 
that the Nelson-Newton law leads uniquely to Schrodinger quantum mechan- 
ics, advocating the view that the Schrodinger equation is a way of linearising 
a very complicated non-linear system of equations equivalent to the Nelson- 
Newton law. 

In this paper we will use Nelson's stochastic mechanics to take the cor- 
respondence limit of the atomic elliptic state for the Coulomb potential and 
show that this limit gives Kepler's laws of planetary motion. 

1.2 Nelson's stochastic mechanics 

We begin with a brief account of Nelson's stochastic mechanics based on the 
summaries in [6, 7]. Let X be the position of a particle of unit mass diffus- 
ing in M.'^ according to an Ito stochastic differential equation with diffusion 
constant e^/2, 

dX{t) = b{X{t),t)dt + edB{t), t > 0, X{0) = x, (1) 

where B = (i?i, i?2, • • • , B^) is a rf-dimensional Brownian motion process on 
a probability space {fi,At,'P), with covariance, 

E{Bi{t)Bj{s)} = 5i,{tAs), 

where E denotes expectation with respect to the measure P. 

We assume that the process X has a smooth density function p{y,t) such 
that for any A E At, 



F{X{t)eA}= I p{y,t)dy. 



'yeA 

Then the density p must satisfy the forward Kolmogorov equation, 

f (2/, t) = V ■ (^ Vp(?/, t) - p{y, t)b{y, t)\ . (2) 



We define the mean forward and backward derivatives D± by, 

f\X{t±K),t±h)-f\X{t),t) 



D±f{X{t),t) :=limE 



X(t) . 



HO [ ±h 

Assuming some mild regularity conditions on /, Ito's formula gives, 

D+/(X(t),t)=(^| + b(X(t),t).V+^A)/(X(t),t). 

In particular, it follows that the mean forward velocity is simply, 

b+{X{t),t) := D+X{t) = b{X{t),t). 
Furthermore, it can be shown that for g = g{X{t),t) and h = h{X{t),t), 

^E {gh} = E {gD_h} + E {hD+g} , 
and it follows, using integration by parts, that for / = f(X{t),t), 

^-/=(|+(&+-eVlnp).V-^A)/. 

This gives the mean backward velocity as, 

b4X{t),t) := D^X{t) = b+{X{t),t) - e^V\np{X{t),t). 

We can also introduce the osmotic velocity u and the current velocity v such 
that, 

u:=^-{b+-b_), «:=^(b+ + b_). (3) 

If we now define the stochastic acceleration by, 

a(X(t), t) = ^ (D+D_ + D^D+) X(t), 
then it follows that, 

a{X{t),t) = (^ + (« ■ V)« - {u . V)u - '^Au^ {X{t),t). (4) 

We now show how this apparently impenetrable expression takes on a deep 
significance in connection with the Schrodinger equation. 



Consider a quantum mechanical particle of unit mass in M^ subject to the 
force field — VV^ where V is some real valued potential. The corresponding 
complex valued wave function ip satisfies the Schrodinger equation, 

in^ = -^A,p + Vij, (5) 

where (letting ip* denote the complex conjugate of ip) the quantum mechan- 
ical particle density is, 

Multiplying through equation (5) by ip* gives the equivalent Schrodinger 
equation, 

ihr^ = -^r^i' + vm\ (6) 

Equating imaginary parts in (6) gives the continuity equation, 

|.v..^o. ,r, 

where j is the probability current. 

If we now write, 

ip = exp(i? + iS), 

where R and S are real valued functions of space and time, then it follows 
that, 

p = exp(2i?), j = hpVS. 

The continuity equation (7) now becomes. 



dt ^ 


-hVSe^^) = V • ( -Ve2« - he' 


Therefore, if we let. 






e^ = h, b = hV{R + S), 


then we have. 


dp ^ A' A 
dt V 2 / 



This is the forward Kolmogorov equation (2) for a diffusion X[t) with forward 
and backward velocities, 

b+{X{t),t)=e^V{R + S), b_{X{t),t) = e^V{S-R). (8) 

A tedious calculation from equations (3), (4) and (8) gives the stochastic 
acceleration in terms of R and S as, 

a(X(t),t) = -V i^-e'— + -{\VR\'- \VS\' + AR)j (X(t),t). (9) 

However, equating real parts of equation (6) gives, 

f = ^(|Vi?P-|V5p + A/2)-^. (10) 

Therefore, combining equations (9) and (10) gives, 

a(X(t),t) = -V\/(X(t),t), (11) 

where V is the potential in the Schrodinger equation. Equation (11) is the 
Nelson-Newton law for a particle of unit mass. 

The Nelson-Newton law argues that to investigate the Bohr correspon- 
dence limit of a quantum state, we should investigate the Nelson diffusion 
process in the appropriate limit. 

1.3 Outline of the paper 

Section 2 starts from Pauli's elegant solution of the hydrogen atom. We 
follow the results of Lena, Delande and Gay [11] to obtain the Schrodinger 
wave function corresponding to an atomic elliptic state in M^ and then find its 
related Nelson diffusion process. In Section 3, we derive the formal limiting 
wave function for this atomic elliptic state according to Bohr's prescription. 
We derive the corresponding limiting Nelson diffusion and present simulations 
of the limiting process showing how the process converges in large times to 
a particular ellipse which we call the Kepler ellipse. We investigate the 
behaviour of the drift for this limiting diffusion and demonstrate that it has 
a finite jump discontinuity across part of the semi-major axis of the Kepler 
ellipse. This singularity is the Bohr correspondence limit of the nodal surfaces 
of the atomic elliptic state wave function. 

For the remainder of the paper we restrict our wave functions and dif- 
fusion processes to the putative plane of motion {z = 0). In Section 4 we 



introduce a new set of elliptic coordinates, the non-orthogonal Keplerian el- 
liptic coordinates, and show how these greatly simplify the limiting Nelson 
diffusion process. 

In Section 5, we analyse the two dimensional limiting diffusion process. 
We discuss the behaviour of the invariant measure for this system and demon- 
strate that it is sharply peaked on the Kepler ellipse. Taking an idea from 
earlier work on excursions in stochastic mechanics [19], we introduce a simi- 
larity transform for the generator of the Nelson diffusion process and use this 
to discuss the convergence over time of the density for the Nelson diffusion 
processes to their invariant measures. We also show that there is a positive 
probability of hitting the drift singularity in a finite time. Thus, to avoid 
inventing an artificial boundary condition, we restrict our path space so that 
our diffusing particle avoids this singularity and we estimate the probability 
for this to happen. On the restricted path space we see that the Bohr corre- 
spondence limit reduces to the underlying deterministic dynamical system. 

In Section 6 we look in detail at this underlying deterministic dynamical 
system. We show that this system has the Kepler ellipse as a stable peri- 
odic orbit and we derive Kepler's laws of planetary motion for a particle on 
this orbit. We identify which paths avoid the singularity and for these paths 
discuss the asymptotic stability of the Kepler ellipse. We show that the in- 
variant density can be used to construct a Lyapunov function for the system. 
Using this, we show that all motions starting outside the Kepler ellipse con- 
verge to Keplerian motion on the Kepler ellipse. For motions starting inside 
the Kepler ellipse the result is more difficult to prove and we leave this to a 
future publication. We conclude by highlighting some surprising symmetries 
within the diffusion process and use these symmetries to show that mild lo- 
cal instabilities occur in the dynamical system for e > 1/a/2 where e is the 
eccentricity of the Kepler ellipse. 

Thus, we see that in M^, for orbits starting outside the Kepler ellipse and 
avoiding the singularity, the long time behaviour of the Bohr correspondence 
limit of Nelson's stochastic mechanics for atomic elliptic states is Keplerian 
motion on the Kepler ellipse. We believe that we can extend this result to 
M^ but so far we have not found a convenient coordinate system in which to 
carry out the analysis. 



2 The atomic elliptic state 

2.1 The wave function for the atomic eUiptic state 

In this section we follow the work of Lena, Delande and Gay [11]. Recall 
that the atomic circular state |Circ)„ corresponding to a Keplerian circular 
orbit is given by, 

(a3|Circ)„ = "^ n,n-l,n-lix) , 

where x = {x, y, z) and '^n,i,m (with / = 0, 1, . . . , n— 1, \m\ < I and m, n G N) 
is the nodal Schrodinger wave function for the Hamiltonian, 

^ |y I 

The vectors p = {pi,P2,P3) and q = (q'l, q'2, ^s) are the momentum and 
position operators in cartesian coordinates for the orbiting quantum particle 
with [qk,Pi] = ihSki for k,l = 1,2,3. We now work in suitable units so that 
h = 1 and fi = 1. Then, for the orbital angular momentum L, where 

L = (Li, L2, L3) = {qxp), L' = {L\ + Ll + L^), 

we have, 

-^ ^n,l,m H^ "•" ^ )^ n,l,rm -'-'3^ n,l,m ''^^ n,l,m) ^ ^ n,l,in 7) 2 "i^i*^' 

The state we consider is, 

|Elliptic(^))„ = exp(-iM2)|Circ)„, (13) 

where, 

A = (A„ A, A,) = -^ gP^L-Lxp) q 

^ ' ' ' ^^2E V 2 \q\ 

is the Hamilton-Lenz-Runge vector on the eigenspace of H with eigenvalue 
E. We will show that on the space where E = En = — l/(2n^), the state 
I Elliptic (^))„ as n I oo corresponds to the elliptic atomic state with eccen- 
tricity e = sin^ for some 6 G (0, vr). 

Pauli proved the fundamental identities, 

[A,, A2] = iLs, [A2, L3] = iAi, [L3, A,] = 1^2, 



showing that Ai,A2, L3 (the constants of motion in the 2-diniensional prob- 
lem) generate the dynamical symmetry group SO (3) for this situation. Now 
define, 

(L3)(^) = (Elliptic(^)|L3|Elliptic(^)), 
{Ai){e) = (Elliptic(^)|Ai|Elliptic(^)). 

It follows that, 

"^^ {L^m = -{L,m, -^{A,m = -{A,m, 



giving, 

(L3)(^) = (n - 1) cos(^), {A,){e) = {n-l) sin(^). 

Moreover, if we consider the classical identity relating the eccentricity e of 
the elliptic orbit to the angular momentum L and the energy E^ 

e^ = 1 + 2L'^E, 

and take the limit as n — > cxd, we find e^ = sin^(^) as asserted. A clever 
argument using coherent state representations and the Kustaanheimo-Stiefel 
transformation yields, for ^0 = arcsin(e) where < e < 1, to within a 
multiplicative constant (after reinstating h and the force constant /i), 

^,,n(a;) := (a;|Elliptic(6'o))„ = exp i- — \x\\ Cn-i{nu), (14) 

where 

/i / X iyVl — e^' 

with X = {x,y,z), En = — /i^/(2A^), e"^ = h, \ = ne^ and Cn-i a Laguerre 
polynomial. 

The nodes of the wave function ip^^n are located on a series of (n — 1) 
hyperbolas in the plane y = each with axis x, eccentricity 1/e and focus at 
the origin (see Figure 1). 

Associated with the atomic elliptic state |Elliptic(^o))n where arcsin^o = 
e, there is a special ellipse which we will call the Kepler ellipse. 

Definition 1. The Kepler ellipse is the ellipse in the plane z = with 
eccentricity e and semi-major axis a with one focus at the origin given in 
cylindrical polar coordinates by, 

a(l — e^) , , 

r=Y^ ^, z = 0, 16 

1 + e cos 

where a = X'^/fx and x = f cos 9, y = fsinO with f = \/x^ + y^. 



^=^(I=-|-T-^^^- ), (15) 




Figure 1: The nodal surfaces in the {x, z) plane for ip^^n for e = 0.5 and n = 5, 
10, 20. 

We will show that the correspondence limit of the Nelson diffusion for 
the atomic elliptic state with wave function ^^ „ is Keplerian motion on the 
Kepler ellipse. 

Remark 1 . Throughout this work the letter e will refer to the eccentricity as- 
sociated with the atomic elliptic state. When we use the exponential function 
we shall always write exp(a;) rather than e^ to avoid confusion. 



2.2 The Nelson diffusion for the atomic eUiptic state 

We now derive the Nelson diffusion process X^'" associated with the wave 
function tp^^n for the atomic elliptic state |Elliptic(^o))n- If we write tp^^n = 
exp(i?e,n + iSe,n) then X"'" satisfies, 



dX^'"(t) = b,,„(X^'"(t)) dt + e dB{t), 



with. 



We will now find the drift 6e,n for the atomic elliptic state. 

The wave function ■?/'e„(a3) satisfies the Schrodinger equation, 

I A . . , . Ai 



--e AV^,,„(x) 



i!e,n{x) = EnAA^)^ 



where e^ = h, En = — /i^/(2A^), A = ne^ and x = {x,y,z). Defining Z^A^) 
by, 

(17) 



Z,Jx) := -ie 



- e V (S'e,„ — iRe,n) 



■^eA^) 



it follows that, 



-y V ■ Z,,„(a^) + ^Zl^ix) -A = E„ 





-1.5 -1 -0.5 0.5 



Figure 2: The x component of the drift b^.n in the plane z = for n = 20 
and e = 0.5 shown as a surface plot and a contour plot (black = — oo, white 

= +oo). 

where in cartesians Z^ „ can be written, 



C-iM\ X ^ /i C'^_i{nu) 



Cn-i{nv) J \x\ \e Cn-iinp) 



(i,-v^ 



e^O,, 



(19) 



z,,„(aj) = y 1 1 - ."~^: : I — + 

and z/ is defined in equation (15). 

Combining equations (17) and (19), the drift for the diffusion X*^'"" is, 

b,,„(aj) = eV(i?,,„ + 5,,„) = ^{Z,^r^{x)) - 53(Z,,„(x)), (20) 

where 3fJ denotes the real part and SJ denotes the imaginary part. 

Clearly, from equation (19), the drift field b^n will have curves of singu- 
larities in the plane y = corresponding to the nodal curves of the wave 
function ifj^^n shown in Figure 1 and also to the point \x\ =0. These singu- 
larities can be seen in Figure 2 which shows the x component of the drift in 
the plane z = 0. The line of singularities corresponds to the intersection of 
the nodal curves with the plane z = 0. As discussed in [3], at each of these 
nodes the drift is infinitely repulsive with the drift components becoming 
plus or minus infinity depending on the direction of approach. This means 
that the nodes will be effectively inaccessible to the diffusion process. As n 
increases the number and density of the nodes will also increase. 

The invariant measure for the wave function ip^^n in the plane 2; = is 
shown in Figure 3 together with simulations of X*^'". The invariant measure 
has a clear peak on an ellipse in the plane z = 0. Away from this plane it 
tends towards zero. 



10 





Figure 3: The invariant measure in the plane z = for the wave function 
ipe,n and simulations of the process X*^'" for n = 20 and e = 0.5. 

3 The limiting atomic elliptic state 

3.1 The limiting wave function 

We now derive the Bohr correspondence limit of the wave function '?/'e,n as 
n — i> oo and e — >■ with A = e^n fixed. Recall the function Z'^ „ from equation 
(17). We define the Bohr correspondence limit of Z'^ „ as, 



-^0,00(33) : = 



lim Z^,n{x), 

nToo,e|0 



where the limit is taken with A as a fixed real number. 

In what follows the Laguerre polynomials £„(x) are defined according to 
the conventions in [5]. 

Lemma 1. Let £„(x) denote the n^^ Laguerre polynomial and X be a fixed 
real number. Then, 



"Too..io Cn-i{nu) 2 \ V 
Proof. For Laguerre polynomials, 

vC'^_^{v) = {n- l)Cn-i{v) - (n - l)Cn-2{v), 

and so, 

C-'n-i{v) n-l n-lCn-2{v) 



jO-n-liv) 



V Cn-l(v) 



11 



Setting V = nu gives, 



lim — -, — - = lim 



However, Laguerre polynomials also satisfy the recurrence relation, 

nCniv) - (2n - 1 - v)Cn~i{v) + (n - l)Cn-2{v) = 0. 

Thus, if the limit, 

Cn-2{nv) 
p = hm 



A— ne'^ 

exists and is non zero, then necessarily p satisfies, 

--(2-z/)+p = 0, 
p 

proving the lemma. D 

Applying Lemma 1 to Z^^n gives, in cartesians. 



(21) 



2A y V uj\x\'z 
where as expected from equation (18), 



2^o,ooW 1^1 - 2A2- 

Proposition 1. The Bohr correspondence limit of the wave function ipe,n for 
the atomic elliptic state gives the formal limiting wave function, 

2A 

71 



M^) = u^{l + \ll--] exp I -^|a;| + ^ ( 1 - ^1 - ^ 



(22) 



Proof. For the vector Zq^^o the function ip^ is the formal wave function sat- 
isfying. 

This wave function corresponds to Bohr's limit but is only an approximate 
solution of the Schrodinger equation. D 

12 



3.2 The limiting Nelson diffusion process 

We now construct the limiting Nelson diffusion X'^ corresponding to the 
limiting wave function ip^. As in equation (17), 

where ipe{x) = exp{R^ + iS'e). The subscripts emphasise the e dependence in 
i?£ and Se arising from the e^. 

Proposition 2. The limiting Nelson diffusion process X'^ corresponding to 
the limiting wave function ip^ satisfies the ltd equation, 

dX'{s) = b{X'{s)) ds + e dB{s), 

where b{x) = {bx,by,bz) in cartesian coordinates with, 

K = ^{{a + (3-l)--{a + f3 + l)^\, (24a) 

2A I e \x\ J 

by = |x{("-/^-l)^^^^-(" + /^+l)^}' (24b) 

b. = -^{a + P + 1)^ (24c) 

and, 



1 {e\x\-x-4X^e/^if + {l-e^)y^ 
" ~ '2V {e\x\ - xy + {1 - e^)y^ ^ 



1 (e|x| -X- 2\^e/fiY + (1 - e^)y^ - 4A^eV^^ 
'2 {e\x\ - a;)2 + (1 - e2)|/2 



(25a) 



S = -2A^ev/r^y 

^ fi{{e\x\-xy + {l-e^)y^)a ^ ' 

Proof. The drift term in the limiting Nelson diffusion process X^ correspond- 
ing to the limiting wave function ip^ is, 

b{x) = e^V{R, + S,) = ?fi{Zo,oo{x)) - S(Zo,oo(a;)). 

A simple calculation gives (for z/ as in equation (15)), 



4 
1 = a + i(3, 

V 

and the result follows. D 
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Figure 4: Simulations of the limiting diffusion X^ with e = 0.5. 



Moreover, we can also find the functions R^ and S^ explicitly, 



R. 



2^ 

2^ 



\n[a^+~^'^] +2\n{{l + af + l3'^) + {I - a)a + 



IJ,\x\ 



S. 



where. 



HTgla + ip) +2arg(l + a + i/3) + -/3(f -a) Pa , 



a 



ii 
A2 



X 



x\ 



p 



IJ,yVT 



e/ A^e 

Our main object of study for the rest of this paper is the limiting Nelson 
diffusion process X*^ in Proposition 2. Figure 4 shows several simulations of 
sample paths for X'^. These paths all converge to the Kepler ellipse. We 
hope to recover Keplerian motion on the Kepler ellipse from the underlying 
deterministic dynamical system X^ found from X'^ in the limit as e ^ 0. 

Remark 2. The process X^ is the formal limit of the process X^'" found by 
letting n — ^ oo whilst e — > with A = ne^ fixed, in the drift term. The 
process X^'" before taking this limit satisfies the Nelson-Newton law. 



1 



{D+D_ + D_D+)X''''{t) = -fi 



X^'"(t) 
X^'"(t)|^ 



3.3 The drift singularity 

Considering equations (25) in Proposition 2, it is clear that a and f3 will 
both have singularities at the point |x| = 0. However, j3 will also have a 
singularity along the surface a = 0. Therefore, the drift could be singular at 
any of these points. 

The former is easy to analyse using spherical polar coordinates. 



X 



r cos 6 sin 0, y = r sin 6 sin 0, z = r cos (f). 
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Then, 

^ = 77 a ■ A\2^f, 2^ • 2zj • 2 A K{e,(j))smesm(j)r 2, 

(^(e — cosy sm0j"^ + (1 — e^j sm t/sm (pj jjra 

as r ~ 0, where C, K are some positive functions independent of r. One 
would normally expect the probability of the diffusion X'^ hitting this point 
singularity in finite time in two or more dimensions to be zero. 

We now consider the behaviour on the surface a = 0. Taking equation 
(25a), we define ai and 0:2 so that. 



a = \l- (v^ + a2)- (26) 

Clearly, a necessary condition for a = is ai — q;| = 0, and so working in 
cylindrical polar coordinates, 

X = f cos 6*, y = f sin 6*, z = z, 

we have, 

2 64AV(l-e)(e + l)f2sin2^ 

ai — a2 = 2' 

/i2 ((e2 cos 2e + e^ + 2) f2 _ 4ef Vr^ + z^ cos ^ + 2e2z2) 

and so o; is only zero on 1/ = 0. Moreover, returning to cartesians. 



a|j/=o 



1 I AX^e/ fi — e\Jx^ + 2;2 + a; 



(4A2e//i — e v^a?2~+^ + x) ' 



(27) 
and so we can conclude, 

a = <^ a; = (x, y, 2;) G S, 

where S is the set. 



-e (4AV/i - 7(16AV/i2 



(x,o,z): ^-^^ <^< Vr^^ 

(2^ 

The set E is shown in Figure 5 together with the Kepler ellipse. This sin- 
gularity is what one would expect if we consider S as the correspondence 
limit of the nodal surfaces of the wave function ■i/'^ „ (see Figure 1). It can be 
shown that this singularity leads to a finite jump discontinuity in the drift. 

15 




Figure 5: The drift singularity (shaded) with the Kepler ellipse for e = 0.1, 
0.5 and 0.9. 



3.4 Restriction to two dimensions 

For the remainder of this paper we shall restrict our wave functions ^^.n and 
ip^ and their related diffusion processes X"'"" and X" to the putative plane 
of motion z = 0. We will consider the full three dimensional problem in a 
future paper. 

Proposition 3. Ifip^^x, y, z) is the three dimensional limiting wave function, 
then its restriction to the plane z = 0, is given by, 



i^e{x,y) 



v^ 1 



2A 



X 



z/ 



/i 



\u 



expl -X^va;' + ^' + ^ 



1- 



V 



where 



/i 



^= ^\ Va;2 + y- 



X 



iy^/l 



Proof. We state without proof that the correspondence limit in three dimen- 
sions satisfies, 



^, 



3— dim I „/,2— dim 

U=0 - We ) 



giving the result. 



D 



Proposition 4. The limiting Nelson diffusion process X^ = (X|,Xp re- 
stricted to the plane z = satisfies, 

dX'{s) = b{X'{s)) ds + e dB{s), 
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Figure 6: Simulations of the 2 dimensional diffusion X^ for e = 0.5. 



where b{x) = {bx,by) with, 



'^ - lxi(«+^-^)^(«+^ + ^)v^ 



by = |;^("-/^-l) 



VT 



(a + 13 + 1] 



Vx'^ + y^ 



(29a) 
(29b) 



and, 



( 



a 



\ 



e^ x^ + y^ — X — AX'^e/fi) + {I — e'^)y'^ 



\ (ey^x^ + y"^ — x) +(1 — 6^)1/2 

I (^e^x^ + y^ -X- 2X^e/i2y + (1 - e^)y^ - 4A V/zi^^ 
(ca/x^ + y"^ — x) +(1 — 6^)1/2 , 



-2\'^eVT^l^l 



P 



Proof. This follows from Proposition 2. 



J(e^WT7^-.)\(l-ev)r^ 



(30a) 
(30b) 

D 



Simulations of the restricted limiting Nelson diffusion process are shown 
in Figure 6. 

4 The Keplerian elliptic coordinate system 

In this section we will define a new two dimensional coordinate system to 
simplify the limiting diffusion X'^. Recall from Definition 1 that the Kepler 



17 



ellipse is the ellipse with eccentricity e and semi-major axis a = A^//x with 
one focus at the origin and the other at (— 2ae, 0). 

In the plane z = 0, the singularity S defined in equation (28) reduces to, 

[a;,0): ^ < a; < I . (31) 

1 + e J 

We want to find a coordinate system which will simplify the complex expres- 
sions a and jS in Proposition 4. From equations (26) and (27) it is apparent 
that on the singularity tti becomes a perfect square. A simple calculation 
shows that the same also happens on the Kepler ellipse. 

In fact if we consider this square root term in polar coordinates, 

8 (2A^ — rX^fi) e^ + rfi cos 6 (8A^ — 2r/i + er/x cos 6) e + r^/x^ 
r2/i2(e cos ^ — 1)^ 

and evaluate this on an arbitrary ellipse of semimajor axis 7 and eccentricity 
c, that is where r = (1 — c^)7/(ccos^ + 1), then we have, 

IGe^c^A^ - 8ec/iA2 (7(0^ - 1) + r{ec + 1)) + fi\e-f{c'^ - 1) + r{c + e))^ 

«! = 2 • 

(e7(c^ — 1) + r(c + e)) /i^ 

The numerator is a quadratic in r and so forms a perfect square when its 
discriminant is zero. This only occurs when 7 = 2ae/(c -|- e). Therefore, 
there is an infinite family of ellipses, including both the Kepler ellipse and 
the singularity (as a degenerate ellipse), on which ai forms a perfect square. 

Definition 2. The family of ellipses Sc are non-confocal ellipses with ec- 
centricity |c|, foci at (0,0) and (— 4aec/(l + e),0) and with semi major axis 
2ae/(e + c) where — e < c < 1. 

Using the ellipses Su we introduce (-u, v) coordinates. 

Definition 3. The Keplerian elliptic coordinates (m, v) in the plane (x, y) 
are defined by, 

2ae(cosf — -u) 2ae^/l — u'^siniv) 

x = . , y = . , 32 

e + u e + u 

where — e < m < 1 and < f < 27r. 

These are non-orthogonal coordinates. The ellipse Sc corresponds to the 
coordinate curve u = c for some constant c G (— e, 1]. The Kepler ellipse is 
the curve u = e (i.e. Sg) and the singularity is the degenerate ellipse given 
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Figure 7: The Keplerian elliptic coordinate system for e = 0.1, 0.5 and 0.9 
with the singularity {u = 1) and Kepler ellipse {u = e) in bold. 

by M = 1 (i.e. Si). The ellipse at infinity is given hj u = —e (i.e. S^e)- The 
curves v = k for some constant k, are similar to hyperbolas. The coordinate 
curves are shown in Figure 7 for several values of e. It is important to note 
the bunching of the curves of constant u which occurs at the left hand end 
of the singularity for large eccentricities. 

The usefulness of Keplerian elliptic coordinates is immediately apparent 
as equations (30) simplify to give, 



a 



v/(i-^')(i 



P 



(e + u) sinw 



1 + e-u — (e + -u) cosf ' ^ 1 + e-u — (e + -u) cosf ' 

and the drift from equations (29) becomes. 



(33) 




eu - 


- {e — u) cosi; ~ {e + u) sint? + ^^{1 — e^) (1 — u'^) — 1 






e(l — ucosv) 


•) 


VI- 


- u^ (e cos V — e sin f + 1) + Vl — e^ {u cos v + u sin v — 


1) 



efl 



■UCOSf 



The singularity S becomes the line u = \ with approaches from above (i.e. 
y > 0) corresponding to < f < vr and approaches from below (i.e. y < 0) as 
TT < f < 27r. In this manner the singularity is opened out onto the boundary 
of our new coordinate space which forms a cylinder. It is also apparent that 
there is a finite jump discontinuity in /9 across S as. 



/5| 



smf 



u=\a]^v 



1 — COS V ' 



/5| 



smf 



u=1,v=2tt—v 



1 — COS V 



We can also use our new coordinates to rewrite the diffusion process X'^ 
in (m, v) space. 
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Proposition 5. The diffusion X^ can he written in terms of the Keplerian 
elliptic cooridinates {u, v) as, 

dX: = MX:,X:){6,(X:,X:)dt-eM(X:,X^)-(d5i, di?2)}, 
where 

N{u,v) = f (e + «)(! — u^) cosf , (e + u)Vl — -u^ sint> j , 
M{u,v) = i (1 + eu) sinf , — Vl — u'^{e + cosf ) 

with, 

h{u,v) — 



2ae{l — ucosv){eu + (e + -u) cosf + 1) ' 
and, 

bu = e^Iu{u,v) - — (e + M)Vl - m^ 

( Vl — e^ (m + cos V — sin f ) — Vl — -u^ (e + cos f — sin f ' 

— — I a/(1 — e^) (1 — u'^)(e + cos v + sin f ) 
^1 + e^) - 2e - (e^ + 2Me + l) cos t; - (l - e^) sin t; j . 
wii^/i i(/ie ltd correction terms, 



by = e Iv{u,v) 



-u[ 



-(e + uf 



^ ^ 4ae{l + eu+{e + u)cosvyy ^ ^^ ^ ^ 

+2m ((e + uY - (1 - m2)(i + eu)) cost; - (1 - m2)(i - e^)], 

:r / X (e + -u)sint; /, ,, , ,, 

4K^) = , .-.^^ ^ / ^ X -^[2iu + ey-il + euf 

+ (e + «)(€« + 1) cosf 1 . 



Proof. Tliis follows from applying Ito's formula to equations (32) which define 
the Keplerian elliptic coordinates. D 
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5 The limiting Nelson diffusion X^ 

5.1 The invariant measure and the infinite time Umit 

We want to show that the (time dependent) density p^{x,t) for the diffusion 
X^ with any given initial distribution Pg(cc) will converge in the infinite time 
limit to a density concentrated on the Kepler ellipse. We begin by showing 
that the invariant measure is this physically correct density in the limit as 

The invariant measure for the Nelson diffusion process X^'" corresponding 
to the atomic elliptic state wave function ■i/'^.n = exp(i?e,n + ^'S'e,„) is given by, 



llV'e.nIP ffexp{2R,^n{x))dxdy' 

This density is shown in Figure 3 where n = 20 and e^ = X/n. For the limiting 
wave function ip^ = exp{R^ + iS^) the invariant measure is pf, where 

.-(^) .= exp{2R,{x)) 

Theorem 1. The invariant density pf coming from the Bohr correspondence 
limit of the atomic elliptic state has the Kepler ellipse £e os a manifold of 
maxima on which it attains a constant global maximum. This maximum is 
sharply peaked on Se as e ^ 0. 

Proof. We can write R^ in terms of Keplerian elliptic coordinates {u, v) giv- 
ing, 

..p^ ..^,'l + e«+v/(l-e2)(l 



e + u 




exp 



2ap lu — e + (e-u — 1 + a/(1 — e^) (1 — u'^) j cosf 



On Se this is constant, 

= 0. 



exp(2i?,(e,t;)) = ( ^ j ^ — exp(2i?,(M,t;)) 



u=e 



Moreover, a simple calculation produces, 

d 

— exp{2R^{u,v)) 
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0, 





Figure 8: p^ in {u,v) and {x,y) space for e = 0.5 and e = 0.1. 

so that £e is a manifold of stationary points. A further calculation shows that 
this is the unique maximum for this function. Clearly, taking into account 
the normalisation factor in p^, the invariant density will become sharply 
peaked as e — ;> 0, as shown in Figure 8. D 

We now want to show that p^{x^t) will converge to p'^{x) in the infinite 
time limit. For this we introduce a similarity transform first developed in 
[19] and applicable to any Nelson diffusion X with diffusion constant e^/2. 
The generator Q of the diffusion X is the operator, 

^ = le2A + 6.V, 

where b = e^V(-R + 5") for suitable real functions R and S such that the 
corresponding wave function can be written ip = exp{R + iS). 

Lemma 2. For the generator Q of a Nelson diffusion process X with diffu- 
sion constant e^/2, 



g = exp(-(i? + S)) -H/e' exp{R + S) 



where H is the formal Hamiltonian, 

and b = e^V(-R + S). Moreover, for ip = exp(_R — S), 
so that ip is the formal ground state for H . 
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Proof. This follows from a calculation using the identity, 
A{fg) = fAg + gAf + 2Vf-Vg, 

together with the real valued nature of E, V, R and S. D 

Recall that the atomic elliptic state wave function ip^^n associated with the 
diffusion X^'" is an exact solution for the Schrodinger equation. We begin by 
showing how for the diffusion X^'" with any given initial distribution p^^^{x) 
the (time dependent) density pe^n{x,t) should converge in the infinite time 
limit to the invariant measure p^^{x). 

The density p^^ni^^t) satisfies the forward Kolmogorov equation (2), 

C'Pe,n[X,t) /->* / ,\ / n\ /' \ 
g-^ = G Pe,n{x, t), PeA^, 0) = Pe,n{^)- 

Therefore, using Lemma 2, we can write, 
pe,n{^,t) = (^exp(i?,,„ + S^A e^P \^-tH,.a/t^j exp(-(i?,,„ + 5',,„))p°„j (aj), 

where H^^n = \{~^^^ + ^en + ^^^ ' ^e,n) IS ^ formal Hamiltonian associated 
with the diffusion X^^n- For fixed n, the work of Blanchard and Golin [3] 
guarantees that with probability one the diffusion X^'" cannot reach the 
nodal surfaces in a finite time because they are infinitely repulsive. Therefore, 
for this case we do not need to worry about specifying a boundary condition 
at these points. 

However, if we are to consider the case n ^ oo which gives the limiting 
diffusion X^, we need to specify a boundary condition for p^{x,t) across the 
limiting singularity S. If we assume that the forward Kolmogorov equation 
is valid on S then, 

-discs ie^—\np^{x,t)j = discs(b-n), 

where n is the unit normal to the singularity S and discs is the discontinuity 
across S. The corresponding boundary condition for 'D{H) (where V denotes 
the domain) reduces to, 

discs (^ A= discs ihy), ^eV{H). 

The putative ground state for H^ = |(— e^A+b^+e^V-b) is ^^ = exp(R^ — S^), 
which satisfies our boundary condition. 
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Assuming that there is a spectral gap, when the hmiting H^ is self adjoint 
and satisfies this boundary condition with a genuine ground state ip^ we 
obtain, 

Pe{x,t) -^ Ceexp(2i?,(a3)), 

as t — !► oo, where 

// p^{x) exp(— 2S'e(x)) dx dy 
""'^ JJeM'2{Re-S,)ix))dxdy- 

This suggests that the Kepler ellipse emerges in the infinite time limit. 

5.2 The limiting diffusion as e ^ 0. 

Using the methods of Veretenikov [20], and Blanchard and Golin [3], we 
can prove the existence and uniqueness of solutions to the Ito stochastic 
differential equations defining the limiting Nelson diffusion X*^ even though 
it has a singular drift field. However, we cannot easily control the limit of 
X^ as e ^ 0. 

We must therefore restrict X'^ to sample paths avoiding the singularity S. 
We do this by estimating the probability that up to a fixed time t the process 
X^(s) avoids the interior of a small ellipse surrounding the singularity. Then, 
using the methods of Freidlin and Wentzell [10], we show how to obtain the 
underlying deterministic system X^ as the limit of X'^ as e ^ 0. 

Definition 4. For small 6 > 0, define the first hitting time, 

r,(int(^i_5)) = inf {s > : X'{s) G int(^i_i), X^(0) = x} , 

where Sis is a small ellipse which surrounds the singularity Si. 

Since mt{Si^s) is open, it follows that, 

P{r,(int(^i_5)) >t} = \imE\-K [ x(int(^i-5))(X^(s)) ds 

kToo [ Jq 

Moreover, for x e ext(£^i_5) the drift b{x) is Lipschitz continuous in space. 
The similarity transform in Lemma 2 leads naturally to the following conjec- 
ture for small e^. 

Conjecture 1. 

P{r,(int(^i_5)) > tj^exp {-{R, + S,)) (exp (-ti/f /e^) exp(i?, + 5,)) (x), 
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where, 

Hf = lim (h, + «:x(int(^i_,)) 

kJoo \ 

is the Dirichlet Hamiltonian with Dirichlet boundary conditions on ellipse 
Si-s, with, 

for b = e^V(-Re + S^), H^ being the self adjoint extension of H^ with domain 

Co~(M2\S). 

For our desired result we want there to be a spectral gap for H^ . In this 
connection it is worth noting that if V ■ b is bounded below and b has a 
unique global minimum at a^min, then, 

^. ~ ^ (-e'A + e^V ■ 6(a;,„in) + b") , 

and according to a celebrated result of Simon [15, 16, 17, 18], a{H^/t^) ~ 
(t{Hq) where, 

1 1 d'^V 1 

2 2 OXiOXj 2 

with V = b'^/2. 

For the limiting diffusion X^, 

fi{e + u) ieu + {e + u) (cos v + sin v) + a/(1 — e^)(l — u'^) + 1 

V3, ■ b = ^^ , 

4aeA(M cos w — 1) (cm + (e + m) cos w + 1) 

(35) 
(where —e<u<l) and, 

/i(e cos v + l) (j{l-e^){l-u^) - l) 

b' = Y( n -■ (36) 

ae"^(-u COS f — 1) 

Moreover, from equation (36), b is symmetric about the y axis {v 1— > —v) 
and so is continuous across the singularity S excluding the point |a;| = 0. It 
can be shown that in polar coordinates, b = 0(r~2) as r ~ uniformly for 
9 G (0, 27r) and also that b has a unique global minimum when u = e^/(2— e^) 
and V = T^ corresponding to the point, 

- ( ~^^ c\ 

^^"'^" V(l + e)(2-e)'7' 
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at which, 



b^fiCmin) = (1 - e) 



2a 

However, from equation (35), it can be shown that V • b is not so well 
behaved. It has a jump discontinuity across S (where u=l) and also blows 
up at each end of S (where v = and v = it). Firstly, in polar coordinates, 

— 1/2 

rV -b^— lasr^O uniformly in 6, but also |V ■ b| ~ r^ as r^ ^ 
but not uniformly in 6 where Vg = \x — (— 4ae/(l + e),0)| , meaning that 
the divergence is not bounded below at the points {x, y) = (— 4ae/(l + e), 0) 
and {x,y) = (0,0). Nevertheless, we hope to pubhsh a proof of the above 
conjecture in the near future using the above and results of Wang [21]. 

Now take a sequence of real numbers tj such that Cj ^ as j ^ oo. 
Then define the sequence of limiting Nelson diffusions X^^ by, 

dX'^{s) = b{X'^{s))ds + ejdB{s), sG(0,t), X'^{0) = x, 

and the process X° by, 

dX\s) = b{X%s)) ds, s G (0, t), X\0) = X, (37) 

where b is as in equations (29). Set, 

Al = {iu : r^(int(^i_5)) > t} , 

where r^ is the first hitting time for X*^^. Then we have the following result 
which is essentially due to Freidlin and Wentzell [10]: 

Theorem 2. If^e'^ < oo, then, 

F{X''{s) -^ X°(s), j -^ oo, uniformly s G (0,t)| A^} = 1. 
Proof. By definition, 

X'^{u)-X^{u)= {b{X'^{s))-b{X\s))) d.s + ejB{u), MG(0,t). 
Jo 

If we restrict u E Al, then for some Lipschitz constant i^ > 0, 

PU 

\X'i{u)-X\u)\<K / \X'^{s)-X\s)\ds + ej\B{u)\. 

Jo 

Set, 

f{t)= f\X^^{u)-X\u)\du, 
Jo 
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then, 

f{u) = \X^^{u) - X\u)\ < Kf{u) + e,\B{u)\, u G (0, t). 
Therefore, 

and so, 

f{s) < ejexp(_ft's) / eic^{—Ku)\B{u)\(\u, 
Jo 

f{s) < ejKexp{Ks) f exp{-Ku)\B{u)\du + ej\B{s)\. 

Jo 

Therefore, 

sup \X'^{u)-X^{u)\<3ej sup \B{u)\. 
se(o,t) se(o,t) 

For any constant c > 0, 

pi sup \X'^ (u) - X\u)\ > c\ < pJ sup \B{u)\>^ 



< 2exp 



C2 



ISe^^t 



Therefore, since exp(— x) < x ^ for x > 0, for any c > 0, 
J^p/sup |X^^H-X°(^)|>cl<^5] 



e^ < oo. 



The result now follows from the Borel-Cantelli lemma. D 

The last result is vacuous unless P {^i} > 0. We can estimate this prob- 
ability with Jensen's inequality giving, 

eLxp(-kJ x(int(£i_5))(X:0ds 

> exp (-kE U x(int(^i_5))(X:0 dsj") . 

If we choose S such that the Lebesgue measure, Leb(int(£^i_a)) = ht{ej)K~^, 
and let K I cxD then formally. 

Pa; {Al} > exp -htitj) / ds sup piix,y) , 



se(o,t) 

yeE 
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where pi is the transition density for X^^ and the supremum is taken over 
y G S^, the upper and lower parts of the singularity, with pi possibly dis- 
continuous across S. 

By methods of Wang [21], we have, 

snpplix, ?/) < - exp ( c I ^ + ^ J J , 

where ej G (0, 1] and t > 0. Then choosing ht{e) = eexp (— c/e^) exp (— ct/e'') 
gives. 



F{Ti>t}+¥{ri<A}< (^j'^ . 



Moreover, the term P{r^ < A} can be made arbitrarily small compared to 
P{r^ > t} for small A, giving the desired result. 

6 The Keplerian dynamical system 

6.1 The dynamical system and Kepler's laws of motion 

In the last section we showed that for paths avoiding the singularity, the lim- 
iting Nelson diffusion X*^ converged to the underlying deterministic system 
X" as e ^ 0. We now consider the behaviour of this deterministic system. 

Definition 5. The Keplerian dynamical system is given by the equations, 

x = b{x) ^ x = h^{x,y), y = by{x,y), 
where b^ and by are as defined in equations (29). 

The vector field for this dynamical system is shown in Figure 9 together 
with the Kepler ellipse and the singularity S. 

Remark 3. Although we use the name 'dynamical system' for the differential 
equations in Definition 5, they do not necessarily give a system satisfying the 
definition of a true dynamical system (see [1, 2] for a full definition). This is 
because the existence of the singularity S means that we cannot use the stan- 
dard methods of differential equations to ascertain the existence, uniqueness 
and extendability of solutions through every point of M^. However, as we will 
show, any path starting outside the Kepler ellipse will stay away from the 
singularity and so we can assume that the definition applied to the exterior 
of the Kepler ellipse will produce a dynamical system. Where necessary we 
will assume that we are working solely on the exterior of S^- We hope to 
extend this work to include all paths using the work of Fillipov [9] . 
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Figure 9: The vector field of tlie Keplerian dynamical system with e 
and 0.9 showing the Kepler ellipse and singularity. 



0.5 



We begin our anlaysis of the Keplerian dynamical system with a very 
simple derivation of Kepler's laws of planetary motion. 

Theorem 3. The Keplerian dynamical system has the Kepler ellipse Se as 
a periodic orbit. Moreover, a particle moving on the periodic orbit Sg obeys 
Kepler's laws of motion with force constant fi and energy E = — /i/2a. 

Proof. Using Proposition 5, we can write the deterministic system in the 
Keplerian elliptic coordinates as, 

u = bu{u,v), v = by{u,v), (38) 

with e = 0. The Kepler ellipse Se has the equation u = e and a simple 
calculation gives 6„(e,f) = so that Se is a periodic orbit. 

Moreover, for u = e, equations (32) which define u and v reduce to give. 



X 



acosv — ae, 



y 



iVl^ 



e^ smf , 



so that V is the eccentric angle for the ellipse Se 
motion (38) reduce to. 



Then, the equations of 



[e,v 



/x 



Aafl 



e cos V 



e cos V 



which is the Kepler law of motion in terms of the eccentric angle v with force 
constant fi and energy E = —fi/2a, where a = X'^/fi. D 

In the following sections we will consider the convergence of trajectories 
of the Keplerian dynamical system to the Kepler ellipse. In the next section 
we will show which initial positions produce trajectories which avoid the sin- 
gularity. In the subsequent sections we will consider the asymptotic stability 
of the Kepler ellipse and we will show that all trajectories which avoid the 
singularity will converge to the Kepler ellipse in such a manner that they will 
also obey Kepler's laws in the infinite time limit. 
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6.2 Avoiding the singularity 

We want to show which initial positions give rise to trajectories that avoid 
the singularity. For the dynamical system, 

x{t) = b{x{t)), t > to, x{to) = xo, 

we denote the solution x{t) = x{t; xo,to). 

From equation (38), the Keplerian dynamical system can be written in 
Keplerian elliptic coordinates as, 

u = bu{u,v), v = K{u,v). 

The drift 6„ tells us about the motion towards and away from the singularity 
E which is given by -u = 1. Thus, we are particularly interested in the sign 
of 6„. 

Lemma 3. For the Keplerian dynamical system represented in Keplerian 
elliptic coordinates, if u E (— e, 1), 

{ u = e, ve (0,27r), or, 

'"-^ ^ \u = J^{v), t;G(7r/2,7r), 

where J^iy) = J^(e,i)(v) / J^{^i^^e)iv) and, 

jF(gj g2)(f) = ei(l — cos f sin f) + 62 (cos f — sinw). (39) 

Proof. A simple calculation from Proposition 5. The curve u = J-'{v) is 
shown in the Cartesian frame in Figure 10. D 

Corollary 1. The singularity S in the cartesian plane is repulsive everywhere 

except for, 

4ae 2ae 
—, <x<-- , 

1+e 1+e 

where for y > the singularity is attractive. 

Proof. For -u = 1, 6„ > for v G (7r/2, it) but &„ < for all other v. D 

Corollary 1 shows that paths which start away from the singularity may 
reach it in a finite time. Using Lemma 3 we can identify which initial positions 
may lead to this behaviour. It can be shown that the curve u = J-'{v) reaches 
its minimum at f = Stt/A and its maximum aXv = 7r/2 and v = n. Therefore, 

1 ir^ ^ ~ 1. ~ -2v^ + 3e 
1 > y-(f ) > e, where e := ^. 

-3 + 2^26 
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e> l/V^"" """" 



Figure 10: The Kepler ellipse £e and singularity (solid line) with the curve 
u = T(y) (dashed line), the ellipse £e (dotted line) and the sign of 6„ indi- 
cated. 

Corollary 2. If Xq G ext(£^e) U ext(£^e) then x{t;Xo,to) will never reach the 
singularity for t > to. 

This tells us that outside an ellipse surrounding the singularity S = £^i, 
the trajectories of the dynamical system will not intersect S. The ellipses Se 
and Sg are shown in Figure 10. 

Remark 4. From this we can conclude the following (see Figure 10): 

1. For e = l/\/2 it follows that e = 1/a/2 so that £e = £e- 

2. For e < l/\/2 it follows that 1 > e > e so that S^ C ext(£^e). 

3. For e > l/\/2 it follows that — e < e < e so that £e C ext(£^e)- 



6.3 The asymptotic stability of Keplerian motion 

We will now look at the attraction of paths which avoid the singularity to 
the Keplerian orbit and also at the stability of this orbit. 

Definition 6. A periodic orbit generating a closed trajectory C is called 
orbitally stable if for any 6i > and any initial position Xq which yields a 
periodic solution traversing C, there exists 62 > such that, 



\xo - Xo\ < 82 



d{x{t;Xo,to),C) < 61, 



for t > to, where d denotes the metric d{x, C) = ini \x — y\ ioi y E C. 
Moreover, if there also exists a ^3 > such that. 



\Xo -Xo\ < 83 



hm d{x{t; Xq, to), C) = 0, 

t— »oo 
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then the orbit is said to be asymptotically orbitally stable. 

For a two dimensional system we have the following result [22] . 

Theorem 4. The closed trajectory C corresponding to a periodic solution 
Xp{t) of period T of a dynamical system is asymptotically orbitally stable if, 




dt< 0. 
Theorem 5. For the Keplerian dynamical system let Xp{t) be a periodic 



trajectory traversing the Kepler ellipse Se with period T = 27r ya^//i. Then, 

db , 




5^^"^ 



^ dt = -27r. 





Proof. We can calculate this integral using Keplerian elliptic coordinates, 

^'^ e cos f + e sin t> + 1 

— ^ dv 

e"' + 2ecosf + 1 

f^ ecosf + 1 

= -2 / dv 

Jq e^ + 2ecosf + 1 

, / _i /I — e fv 

= — lim [ V + 2 tan tan - 

vU\ \l + e \2 

= -27r 
as 1 > e > 0. n 

Corollary 3. The Kepler ellipse Se is an asymptotically orbitally stable pe- 
riodic orbit for the Keplerian dynamical system. 

We can also demonstrate that the Kepler ellipse is asymptotically or- 
bitally stable by using the invariant density p^ from equation (34) as a 
Lyapunov function. 

Theorem 6. Let VlpvI^^) be a real valued function defined in an open neigh- 
bourhood N{C) of a compact set C . Assume that, 

1- ^Lpv(3^) i'S continuously differentiable, 

2. VLpv(a3) is positive definite in N{C) \ C , 
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Figure 11: The Lyapunov function Vlpv in {u,v) and {x,y) coordinates for 
e = 0.5. 

3- VVlpv(3^) ■ b{x) is negative definite in N{C) \ C, 

4. Vlp,{x) = VVlpv(3J) • b{x) = OforxeC. 
Then the compact set C is asymptotically stable. 

The function Vlpv in Theorem 6 is called a Lyapunov function. 
Theorem 7. The function, 

Vipv= {l6/e^)^-exp{2R,)l^,, 

satisfies the conditions of Theorem 6 for the compact set C = Se and the 
neighbourhood N{C) = mt{£_e+Si) \ ^^^{^1-52) for any 61,62 G (0, 1 + e) such 
that 61 + 62 <l + e. 

Proof. This follows simply from Theorem 1. The function Vlpv is shown 
in Figure 11. It is continuous but is not smooth across S. The derivative 
V^Lpv • b is shown in Figure 12. We exclude the ellipse at infinity £^_e, as the 
derivative is zero in the limit. D 

Asymptotic orbital stability tells us that not only does any trajectory 
sufficiently close to Sg converge in the large time limit to Sg but that in this 
large time limit any converging trajectory will move exactly as a trajectory 
on Se- Therefore, the motion of a particle which converges towards the ellipse 
Se must converge to Keplerian motion on Se- 
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Figure 12: The derivative Va;VLpv ■ b{x) in {u,v) and {x,y) coordinates for 
e = 0.5. 

6.4 Region of attraction to the Kepler ellipse 

We now look at which initial positions produce trajectories which are at- 
tracted to the Kepler ellipse in the infinite time limit. As in Section 6.2 we 
are particularly interested in the sign of 6„ as the Kepler ellipse is given by 
the coordinate curve u = e and we know that &„ = on this curve. 

Theorem 8. Let < e < 1/a/2- For any initial position Xq G ext(£^e) o-nd 
any given 6 > there exists a finite time T > such that, 



d{x{t;Xo,to),£e) < 5, 



for all times t > T. 



Proof. For the curve u = J-'{v) to cut the Kepler ellipse m = e we would 
require e < e or equivalently e < l/\/2. Therefore, if < e < l/-\/2 then (as 
shown in Figure 10): 

1. for —e<u<e, hu > 0, 

2. for e < M < J^{v), hu < 0. 

Since b^ > for all u G (— e, 1) and v G [0, 27r) the result follows. 



D 



From the proof of Theorem 8 we see that for < e < l/-\/2, the Kepler 
ellipse is a stable orbit in the (u,f) frame as the curve u = J-'{v) does not 
intersect the Kepler ellipse. However, when e = 1/v^ the curve u = J-'{v) 
touches the Kepler ellipse and for e > 1/v^ it intersects the Kepler ellipse. 
Thus, there is a portion of the Kepler ellipse which is unstable in this frame. 
Despite this we cannot conclude that the Kepler ellipse will be unstable 
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in the cartesian frame as the properties of stabihty and instabihty are not 
necessarily coordinate invariant. As we showed in Corollary 3, the Kepler 
ellipse is actually stable. As can be seen in Figure 13, simulations of the 
dynamical system starting outside the Kepler ellipse always converge to the 
ellipse. The region bounded by the curve u = J-'{v) is where this convergence 
is particularly slow. It would seem sensible to investigate whether e = l/\/2 
is a critical value of the eccentricity. It could be that this apparent instability 
is a property of the {u, v) coordinate system. As we shall see in the next 
section this is not the case! 

If we assume e < l/v2 then we can avoid these problems and combining 
our results we can conclude: 

Theorem 9. For < e < 4^, any orbit of our Keplerian dynamical system 
with a start point outside the ellipse Se settles down in the infinite time limit 
to Keplerian motion on the Kepler ellipse Se- 

Proof. This follows from asymptotic stability combined with Theorems 3 and 
8. The region of convergence is shown in Figure 14 (a). D 

We now consider ways to extend this domain of attraction to include 
eccentricities e > l/v2- Arguing exactly as above we have the following 
result which restricts the location of the trajectory after a finite time to a 
large neighbourhood of the Kepler ellipse. 

Theorem 10. Let 1 > e > 1/a/2. For any initial position Xq E ext{Se) there 
exists a finite time T > such that, 

x{t;Xo,to) e int(£^e) next(£^e), 

for all times t > T. 

Proof. Follows as for Theorem 8 but now e > e so that Se C Sg. This region 
is shown in Figure 14(b). D 

This still leaves us to prove that trajectories in this elliptic annulus con- 
verge to the Kepler ellipse. This can be done using the Lyapunov function 
in Theorem 7 and results of [1]. 

6.5 Surprising symmetries and mild instabilities 

The appearance of e = 1/V^ as an important eccentricity could easily be dis- 
missed as a property of the Keplerian elliptic coordinate system. Numerical 
experiments confirm that although 6„ does switch sign in a neighbourhood 
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Figure 13: Simulations of the Keplerian dynamical system with e = 0.1, 0.5 
and 0.9. 




£-e. 




Figure 14: Theorems 9 and 10: (a) The region of attraction (shaded) for £(. 
with e < 1/v^- (b) For e > 1/v^ any trajectory starting outside the Kepler 
ellipse must converge to somewhere in the shaded region. 
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of the Kepler ellipse for e > \I^J2^ the Cartesian dynamical system still al- 
ways moves towards the Kepler ellipse in this region. This discrepancy can 
be explained easily by considering the bunching of level surfaces of u which 
occurs in the Keplerian elliptic coordinates (see Figure 7) and also the slow 
rate of convergence in this region (see Figure 13), particularly for large ec- 
centricities. It is therefore surprising to encounter this critical eccentricity 
through other calculations. 
Recall that if -u 7^ e then, 

6„ = ^ u= /"^''^^ ,\ , t;G(7r/2,7r), 
^{-\,-e)[v) 

where T{ex,e-2.){y) is as defined in equation (39). 

If we consider the divergence of the drift field calculated in Cartesians 
but presented in Keplerian elliptic coordinates for simplicity (see equation 
(35)) then we find, 

V-6 = ^ «=,iiliz!^, t;G(7r,37r/2). 

Therefore, the curve V-b = is the reflection in the y axis of the curve &„ = 0. 
Moreover, this leads to the conclusion that for e < l/v2, the curve V ■ h 
does not cut the Kepler ellipse, and so on all parts of the ellipse V ■ b < 0. 
However, for e > 1/v^ there is a portion of the ellipse where V ■ b > 
suggesting a local mild instability in this region. 

These instabilities are not just properties of the two dimensional restric- 
tion. Consider the full three dimensional system for small values of \z\. The 
drift in the z direction is given by, 

6. = -|^(a + /3+l)^. 

Thus, ifa + Z^ + l > 0, then any trajectory will be attracted to the 2; = 
plane where it will remain stable. Let us consider the behaviour of a + /9 + 1 
in the plane z = 0. Working from equations (33) it is simple to show, 

a + /5+l = ^ u= },^'~^^,\ , t;e(0,7r/2), 

^(-i,e)(-^^) 

which is a curve in the same family again. It again meets the Kepler ellipse 
when e = 1/a/2. Therefore, when e > 1/a/2, there will be a region of the 
Kepler ellipse where the orbit is unstable in the z direction as within this 
curve q; + /9 + 1 < which will mean 62 > for small z > and 62 < for 
small 2; < 0. These three curves are shown together in Figure 15. 
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Figure 15: The curves 6„ = (long dash), V-b = (dotted) and a+P + 1 = 
(dash-dot hne) with the Kepler ellipse and singularity for e = 0.5 and 0.8. 

In fact the curve a + (3 + l = 0isd. slice of a three dimensional surface 
which, together with the singularity, bounds a region of space where bz will 
be directed away from the plane z = 0. This creates a general instability in 
orbits for e > 1/a/2 as this region is cut by the Kepler ellipse. The effect 
of this in the deterministic case is shown in Figure 16. There is a clear blip 
in the value of z where the particle moves away from the plane z = 0. This 
blip is repeated each time the particle passes through this region, but each 
time the blip becomes smaller as the particle moves closer to the plane z = 
overall. In the infinite time limit the particle still appears to converge to 
the plane z = 0. However, in the stochastic case (e ^ 0) the noise creates a 
displacement such that the particle never appears to settle down into a small 
neighbourhood of the plane z = (see Figure 17). This is caused by the 
presence of the noise, but is apparently exacerbated by these repeated blips 
in bz, which unlike the deterministic case do not decrease in magnitude in 
the infinite time limit (see Figure 18). This effect becomes more pronounced 
as e — >■ 1 since the trajectories pass closer to the origin where b ~ ±r~^/^. 
The periodic spikes in bz shown in Figure 18 are the result of this effect and 
cause the particle to be strongly forced into the plane z = when at the 
perihelion. The blip in bz occurs immediatly after the particle has passed the 
perihilion forcing it away from the plane z = again. 

7 Conclusions 

We have shown that in the Bohr correspondence limit in two dimensions 
the Nelson diffusion process corresponding to the atomic elliptic state for 
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Figure 16: Simulation for e = and e = 0.9 highlighting the periodic blips 



m z. 




Figure 17: Simulation for e = 0.05 and e = 0.99 with the invariant measure. 




Figure 18: The z value (thick) with velocity hz (thin) and driving noise 
S3(dashed) for e = 0.05 and e = 0.99 with a close up of a blip. 
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the Coulomb problem reduces to Keplerian motion on the Kepler ellipse for 
motions starting outside this ellipse. This solves a long standing problem in 
quantum mechanics. In the quantum mechanical setting, Kepler's laws of 
planetary motion need to be augmented with a caveat about the mild local 
instabilities appearing on the elliptical orbit for eccentricities greater than 
l/\/2 which may be experimentally detectable. In any case these results 
merit further study in the setting of planetesimal diffusions. 

We are currently investigating the three dimensional problem which ap- 
pears to be far more difficult to analyse in the absence of any convenient 
coordinate system. 
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